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ABSTRACT 



Aims. In the last four decades it has been observed that solar flares show quasi-periodic pulsations (QPPs) from the lowest, i.e. radio, to 
the highest, i.e. gamma-ray, part of the electromagnetic spectrum. To this day, it is still unclear which mechanism creates such QPPs. 
In this paper, we analyze four bright solar flares which show compelling signatures of quasi-periodic behavior and were observed 
with the Gamma-Ray Burst Monitor (GBM) onboard the Fermi satellite. Because GBM covers over 3 decades in energy (8 keV to 
40 MeV) it can be a key instrument to understand the physical processes which drive solar flares. 

Metliods. We tested for periodicity in the time series of the solar flares observed by GBM by applying a classical periodogram 
analysis. However, contrary to previous authors, we did not detrend the raw light curve before creating the power spectral density 
spectrum (PSD). To assess the significance of the frequencies we made use of a method which is commonly applied for X-ray binaries 
and Seyfert galaxies. This technique takes into account the underlying continuum of the PSD which for all of these sources has a 
P(f) ~ /"" dependence and is typically labeled red-noise. 

Results. We checked the reliability of this technique by applying it to a solar flare which was observed by the Reuven Ramaty High- 
Energy Solar Spectroscopic Imager (RHESSl) which contains, besides any potential periodicity from the Sun, a 4 s rotational period 
due to the rotation of the spacecraft around its axis. While we do not find an intrinsic solar quasi-periodic pulsation we do reproduce 
the instrumental periodicity. Moreover, with the method adopted here, we do not detect significant QPPs in the four bright solar flares 
observed by GBM. We stress that for the purpose of such kind of analyses it is of uttermost importance to appropriately account for 
the red-noise component in the PSD of these astrophysical sources. 

Key words. Sun: flares, Methods: data analysis. Methods: statistical. Methods: observational 



1. Introduction 

. Over the last 40 years quasi-periodic pulsations (QPP) in 
■ solar flares have been reported from observations across 
the electromagnetic spectrum, i.e. from radio waves up to 
the high energetic gamma-rays ranging from sub-second 
timescales up to several minutes (e.g . Parks & Winckler 1969: 
jOfman & Sui! 120061: iLi & Gad 120081: Eakariakov & Melniko^ 
l2009i:iNakariakov et alj2010 a). While there seems to be an over- 
whelming amount of observational data, the underlying physi- 
cal mechanism which could generate such QPPs remains still a 
mystery. Possible processes that are considered are modulation 
of ele ctron dynamics by magneto hydrodynamic (MHD) oscilla- 
tions dZaitsev & S tepanov 1982^, periodic triggering of energy 
releases by MHD waves (Foullon et al. 2005; Nakariakov et al. 
l200d ). MHD flow overstabilities (lOfman & Sui 2006) and oscil- 
latory regimes of magnetic reconnection (IKliem et al.,.20 00). 

In this paper we will present time series and periodogram 
analyses of four solar flares which show an intriguing quasi- 
periodic behavior in their light curves. All of these solar flares 



were observed by the Fermi Gamma-Ray Burst Monitor (GBM, 
iMeegan et al. 2009). GBM is one of the instruments onboard 
the Fermi Gamma-Ray Space Telescope dAtwood et al.l l20()9l) 
launched on June 11, 2008. Specifically designed for gamma- 
ray burst (GRB) studies, GBM observes the whole unocculted 
sky with a total of 12 thallium-activated sodium iodide (Nal(Tl)) 
scintillation detectors covering the energy range from 8 keV 
to 1 MeV and two bismuth germanate scintillation detectors 
(BGO) sensitive to energies between 150 keV and 40 MeV 
(Meegan et al. 2009). Therefore, GBM offers superb capabilities 
for the analyses of not only GRBs but solar flares as well. 

The analysis and interpretation of power spectral density 
(PSD) of solar flares is, in general, difficult. A variety of as- 
trop hysical sources (such as X-ray binaries, Seyf ert galaxies 



^ 



— . Lawrence et al l 119871 iMarkowitz et all l2003h and GRBs 
lUkwatta et all 120091: ICenko et al.l 120101) '). show erratic, aperi- 
odic brightness changes. Solar flares exhibit similar aperiodic 
variations with the general time profile being a sharp impulsive 
phase followed by a slower decay phase. Solar flares, together 
with many other astrophysical sources, thus have very steep 
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power spectra in the low-freque ncy region. This type of variabil- 
ity is known as red-noise (e.g. lGrothl[l975t iDeeter & Bovntorj 
[1982; Vaughan 2005). When determining the significance of 
possible periodicities in the PSD the red-noise has to be ac- 
counted for in order not to severely overestimate t he significance 
of identified frequencies dLachowicz et al.ir2009h . In this paper 
we will account for the red-noise properties when performing a 
periodogram analysis. 

This paper is organized as follows. In Sect|2] we briefly 
present the methodology of the time-series analysis. We provide 
an overview of the red-noise properties in astrophysical sources 
and demonstrate the importance of the red-noise when estimat- 
ing significances. In Sect[3]we will present light curves and peri- 
odograms of solar flares which were observed by GBM. Finally, 
in Sect|4] we will summarize and conclude. 



2. Analysis of data governed by red-noise 



As it was already pointed out by [Mandelbrot & WallisI (1 19691) 
and lPressI ( Il978l) . the human eye has a tendency to identify pe- 
riodicities from purely random time series, i.e. where sinusoidal 
variations are not statistically real. According to Press (197^ 
the strongest eye-apparent period in (actually non-periodic) data 
will be about one-third the length of the data sample. According 
to these authors "three-cycle" quasiperiods should be taken with 
a grain of salt. 

Solar flares fall into the group of astrophysical sources 
where red-noise is important. Red-noise has nothing to do with 
measurement errors or systematics of the detectors, which are 
also called noise. Red-noise is an intrinsic property of the ob- 
served source and is due to erratic, aperiodic brightness changes. 
Contrary to white noise, which displays a flat spectrum in a PSD, 
i.e. is power independent of frequency, red-noise is character- 
ized by a power law of the form of P = Nf'". As a first order 
approximation, red-noise is the realization of a linear stochastic 
and weakly non-stationary process. This red-noise component 
makes the interpretation of the significance of a peak in the PSD 
more complex. 

One way to estimate the significance of induced frequencies 
on top of an underl ying re d-noise continuum in a PSD, was pre- 
sented by Vaughanr (l2005l) . 

In short. iVaughanI suggests to calculate the periodogram nor- 
malized so_fliatthejinits_of_20W^^ (rms/mean)^Hz ' (e.g. 
ISchustedlT89l iPress & RvbTckT 1989). Then, the periodogram 
is converted to log-space in both frequency and power For such 
a log -periodogram one can than clearly identify the power-law 
component in the low-frequency range and a "cutoff"' where 
white noise or an a dditional noise component takes over (e.g. 
lUkwatta et al.ll2009h . One can easily determine the power-law 
parameters by fitting a hnear function to the low-frequency pe- 
riodogram bi ns using the l east-squares method. In this paper the 
method of Vaughan (2005) was slightly modified in that we use a 
broken power law (BPL) to fit the PSD instead of a single power 
law. 



2.1. Red-noise simulation 



It is common practice (llnglis et al.l 120081: iNakariakov et al.l 
1201 Oah to suppress the low-frequency component by de-trending 
the light curves of solar flares. This can be achieved by smooth- 
ing the light curve with a moving average or by applying a 
Gaussian filter and perform the periodogram analysis on the 
residual emission, i.e. the smoothed version is subtracted from 
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Fig. 1. Upper panel: Synthetic red-noise time series with a 
power law index of a - -1.8. Simple (boxcar) moving aver- 
age is overplotted (dashed line). The inset presents the residual 
signal. Middle panel: Lomb-Scargle Periodogram of the resid- 
ual light curve. The 3cr CL has been denoted by a dashed line. 
Lower panel: PSD with the best fit power law (solid green line) 
and the 3cr significance level (dotted line). 



the original data-set. This can give rise to misleading results as 
we will show in the following. 

By randomizing the phase as well as the amplitudes, 
iTimmer & Koenia (Il995h introduced an algorithm to generate 
a purely random time series which shows a Ijf" dependence 
in the PSD. We created a time series consisting of 300 data 
points, evenly spaced by 1 s with a /"' '^ periodogram shape. 
From this light curve, we subtracted a simple moving aver- 
age of 50 s (see Fig. [U. We then calculated the PSD (first de- 
scribed by Lomb (1976 ) and Scargle ( 1982)) and then refined by 
IPress & Rvbickiiri989h on the residuals. The result is striking. 
Although we started with a purely random, red-noise dominated 
time-series we obtain a PSD with three frequencies whose power 
exceeds the 3o- confidence limit (calculated according to Scargl^ 
1982). This approach in signal processing clearly returns false- 
positive frequencies with periods of f «; 33 s, P a; 55 s and 
P X! 100 s, respectively. When the periodogram is calculated 
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over the original light curve, using the l VaughanI (l2005h approach, 
i.e. with a single power law fit to the data, the spectral peaks re- 
main below the 3cr threshold. 

Therefore, by this example, we strongly discourage subtract- 
ing smoothed versions of raw light curves when looking for in- 
trinsic frequencies in the red-noise dominated PSD. 

2.2. Detour: QPPs in GRBs 

The same procedure of light curve reprocessing with a sub- 
sequent PSD analysis has also been applied for GRBs. These 
events are the most luminous flashes of y-rays known to 
mankind, believed to originate from highly relativistic out- 
flows from a compact source wi th Lorentz fa ctors T > 100. 
In 2009, the Swift-M r satellite dGehrels et al. 2004) observed 
ORB 090709A (.Morris et all 120091) . Soon after the detection, 
iMarkwardt et al.l (l2009t) claimed to have found a very un- 
usual behavior, not observed in any other GRB so far, namely 
a QPP like behavior with a periodicity of P w 8.1 s at 
the 12 cr level of significance. This QPP was subsequently 
confirmed and found to be in phase in the data of the 
Anticoincidence System (A CS) of the spect rometer SPI on board 
the INTEGRAL satellite (Gotz et al. 2009), the Konus-WIND 
instrument ( Golenetskii et al. 2009) and the Suzaku Wide-band 
All-sky Monitor (WAM) (lOhno et al.ll2009l) . These latter insti'u- 
ments operate in the energy ranges 80 keV-10 MeV, 20 keV- 
1150 keV and 50 keV-5 MeV, respectively. Swift-BAT, on the 
other hand, is sensitive in the energy range 15 keV-150 keV. 

However, soon thereafter, I Cenko et al.l (1201 Oh showed that 
the interpretation of this QPP strongly depends on the assump- 
tion of the underlying continuum. If, in fact, it is accounted for, 
the significance of the claimed periodicity drops below the 3 cr 
confidence hmit. T his analysis was indepe ndently repeated by 
llwakiri et al.l (1201 Ol) and lde Luca et al J (1201 0) who also took into 
account the red-noise component in the PSD, and only find a 
marginally significant periodicity at the 3cr confidence limit. 

In conclusion to this detour, we emphasize once more the 
importance to account for the red-noise component in the PSD. 
Additionally, we draw attention to the fact that a potential quasi- 
periodic signal is not necessarily significant even if it is identi- 
fied in several instruments with different (but overlapping) en- 
ergy ranges and observed to be in phase across these bands. 



2.3. Method testing 

The Reuv en Ramaty High -Energy Solar Spectroscopic Imager 
(RHESSI, iLin et al.l |2002|) rotates around its spin axis which 
is always pointed towards the Sun. The period of this rotation 
is P ~ 4 s and a PSD analysis of the RHESSI light curves 
is ex pected to display this instrumental signal. We applied the 
IVau ghan (2005^ test to a solar flare which has been observed 
by RHES SI on January 1st. 2005 and where QPPs have been 
reported dNakariakov et al.l 1201 Oal) . This solar flare peaked at 
00:31 UT at a GOES level XI. 7, from the NOAA active re- 
gion 10715 located on disk at N03E47. We used RHESSI data in 
the energy range from 50 keV-100 keV with a fine time resolu- 
tion of 0. 1 s (see upper panel of Fig|2]) and performed two peri- 
odogram analyses on this light curve in the range between 1660 s 
and 1820 s. The first periodogram an alysis was per form ed using 
the cla ssical approach introduced by Lomb| (1 19761) and Scargle 
(Il982|). Similarly to what has b een commonly done in the past 
(e.g. Ilnglis & Nakariakovll2009l) a periodogram has been calcu- 
lated on the residual emission after a simple moving average of 
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Fig. 2. Upper panel: Summed and background subtracted light 
curve of the solar flare observed on January 1st, 2005 by 
RHESSI. The inset shows the residual emission after subtracting 
a simple (boxcar) moving average. Middle panel: Periodogram 
analysis performed on the residual emission. Lower panel: PSD 
with best fit broken power-law (solid green line) and the 3cr sig- 
nificance level (dotted line) indicated. The arrow points to the 
significant frequency at / » 0.244 Hz. 



50 s has been subtracted from the raw data (see middle panel 
of Fig|2]i. With this method, several peaks are found above the 
3 cr threshold. The peak with the highest value of normalized 
power is located at f ~ 0.025 Hz, corresp onding to the period- 
icity reported bv lNakariakov et al] (1201 Oal) . Another peak which 
is worth mentioning is located at / ~ 0.24 Hz which is the ex- 
pected rotational frequency of RHESSI around its spin axis. 

As we will show in the following, the significance of the peak 
at / ~ 0.025 Hz is highly overestimated by the latter method. 
We show a PSD which was calculated using the raw and und e- 
trended light curve applying the technique by Vaughan (/005) in 
the bottom panel of Fig|2] Analogically to the Lomb and Scargfl 
periodogram analysis, we found a significant spectral feature at 
w 0.24 Hz which is the expected rotation period of the RHESSI 
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Fig. 3. Upper panel: Light curve of the giant flare of SGR 1806- 
20 observed by SPI-ACS. Lower panel: PSD of the SGR 1806- 
20 light curve. Best fit broken power-law to the PSD continuum 
(solid green line) and the 3cr significance level (dotted line) are 
indicated. The arrow indicates the recovered rotation period of 
P ~ 7.56 s. Other peaks above the 3cr levels are harmonics. 



instrument. However, this PSD is lacking any other frequency 
above the 3 cr confidence limit. The discrepancy between the two 
methods is easily explained. Firstly, the whole raw light curve 
was used without artificia lly detre nding it beforehand. Secondly, 
the method bv .Lomb. and lScargle, assumes a white noise contin- 
uum and does not take into consideration the red noise compo- 
nent of the solar flare. However, the latter is taken into account 
by the method of Vaughan. In conclusion, we could not confirm 
the reported QPP in this solar flare. We are confident that this 
latter method can be used for further analysis since we believe 
that it is more adequate for sources which are dominated by red- 
noise and the inherent rotational frequency of the instrument is 
found to be significant. 

An additional check was performed using SPI-ACS data of 
the outburst of the soft gamma repeater SGR 1806-20 observed 
on Dec ember 27th, 2004 in th e energy range from 80 keV to 
8 MeV dMereghetti et al.ll2005l) . SGR 1806-20 is known to have 
a rotational period of 7.56 s and the methodology adopted here 
should be able to recover this periodicity. We removed the very 
bright initial pulse and focused on the emission from « 50 s to 
175 s (see Figl3]l. We unambiguously recover the main pulsation 
period (P ^ 7.56 s) together with the first, second, third and fifth 
harmonic. 



We conclude that the here applied methodology is appropri- 
ate and reliable for the further analysis. 

3. Solar flares observed by GBM 

3.1. Solar flare on February 24, 201 1 at 07:29:20. 71 UT 

For the purpose of its analysis we use CSPEC and CTIME data 
of detectors Nal 3, Nal 4 and Nal 5 with a time-resolution of 
1.024 s (4.096 s pre-trigger) and 0.064 s (0.256 s pre-trigger), 
respectively. In the energy range from 50 keV to 1 MeV this so- 
lar flare lasted for about 500 s. The light curve (Fig|4]i consists of 
several peaks and a compellingly looking quasi-periodic behav- 
ior lasting to x 500 s. After de-trending the raw light curve with 
a simple moving average (50 s) the QPP pattern becomes more 
visible ( see inset of FiglUl. Appl ying a standard periodogram 
analysis (iLombl 1 976l : IScarglel 1 98 2) on the detrended light curve 
several peaks are above the 3cr confidence limit as can be seen 
in the middle panel of Fig|4] 

However, a pplying the previously introduced method 
(IVaughanll2005l) on the original light curve, thus taking into ac- 
count the red-noise properties of the source we did not find any 
significant QPP during the solar flare life-time, as is pointed out 
in the lower panel of Fig|4] We performed this analysis making 
use of the CTIME data with a time resolution of 64 ms. The PSD 
has been calculated for the signal spanning from (To - 10) s to 
(To + 600 s), where To denotes the time of trigger. 

3.2. Solar flare on June 12, 2010 at 00:55:05 UT 

For the analysis we use again CSPEC and CTIME data of de- 
tectors Nal through Nal 5. The count rate is not increasing 
significantly for the first 30 s in the 50 keV to 1 MeV energy 
range. After this time there is a sharp increase in the flare bright- 
ness which then decays again very rapidly after 60 s. The whole 
duration of the solar flare in this energy range is approximately 
120 s. Overlaid on top of the observed light curve, again one 
can identify a compellingly looking QPP behavior with a pe- 
riod of ~ 15 s (see upper panel of Fig|5]l. This periodicity is 
apparently significant when applying a standard Lomb-Scargle 
periodogram (see middle panel of Fig|5]l on the detrended light 
curve. The smoothing length was 15 s. 

However, the PSD, computed for CTIME data with a time 
resolution of 64 ms in the interva l (To - 25) s to (Tq + 150) s 
and according to iVaughanI (l2005h . does not show a significant 
periodicity. 

3.3. Solar flare on March 15, 201 1 at 00:21:15.69 UT 

For the analysis we use again CTIME data of detectors Nal 
through Nal 5 in an energy range covering 50 ke V to 1 MeV with 
a time resolution of 0.256 s. The total duration of the solar flare 
in this energy range is approximately 50 s (see Fig|6]l. For the 
Lomb-Scargle periodogram analysis (see middle panel of Fig|6ll 
the light curve was detrended with a simple moving average of 
5 s. Repeating the PSD analysis on the original data set between 
(To - 30) s and (Tq + 100) s according to .VaughanC2005.) we did 
not find any significant QPP. 

3.4. Solar flare on March 14, 2011 at 19:50:17.3 UT 

For the analysis we use the same data type, energy range and 
time resolution as for solar flare #4 in Sect l3.3l The total duration 
of the solar flare is approximately 150 s (see Fig|7]). Contrary to 
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Fig. 4. Upper panel: Summed and background subtracted light 
curve of the solar flare observed on February 24, 201 1 by GBM. 
The inset shows the residual emission after subtracting a sim- 
ple (boxcar) moving average. Middle panel: Periodogram anal- 
ysis performed on the residual emission. Two peaks are above 
the 3cr confidence limit (dashed line). Lower panel: PSD of the 
same solar flare. Best fit broken power-law to the PSD contin- 
uum (solid green fine) and the 3cr significance level (dotted line) 
are indicated. 



the a Lomb-Scargle periodogram which finds a significant peri- 
odicity in the detrended light curve (smoothing length of 10 s), 
the PSD, determined in the time inter val betwe e n (Tn - 25) s and 
(To + 150) s and done according to IVaughanI (l2005h . does not 
find any significant QPP 
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Fig. 5. Upper panel: Same as Fig|4]for the solar flare observed 
on June 12, 2010. Middle panel: Same as FigE] Lower panel: 
Same as Fig|4] 



4. Summary & Conclusions 

We have analyzed light curves of five solar flares observed by 
both GBM and RHESSI. We tested the data for the presence and 
significance of QPPs accounting for the overall shape of the PSD 
by applying the method introduced bv lVaughanI (l2005h . First of 
all, this technique was validated and tested by applying it on raw 
data of the RHESSI satelHte. Any RHESSI fight curve has an in- 
herent period caused by the rotation of the space craft around its 
axis. With the method adopted here we successfully retrieve this 
well known 4 s period. However, we were not able to confirm 
the previously reported QPP of 40 s in the very same solar flare 
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Fig. 6. Upper panel: Same as Fig|4]for the solar flare observed Fig. 7. Upper panel: Same as Fig|4]for the solar flare observed 
on March 15, 2011. on March 14, 2011. Middle panel: Same as Fig|4] Lower panel: 

Middle panel: Same as FigH] Lower panel: Same as FiglH Same as FiglH 



JNakariakov et al.ll2010 ah. An additional check was performed 
by applying the method to SPI-ACS data of the giant flare of the 
well known SGR 1806-20. This magnetar has a rotation period 
of 7.56 s which, together with several harmonics, could be recov- 
ered unambiguously with the procedure presented here. These 
two tests gave us confidence that the method is appropriate to 
test the significance of QPPs in red-noise dominated solar flare 
time series. 

The routine was then applied to four solar flares ob- 
served by GBM. Although all of these solar flares showed 
very promising quasi-periodic features in their detrended light 



were 
ir 



si gnificant (> 3 cr). 



curves at first none of these 

Previous authors (e.g. Inglis et al. 2008; 'inglis & Nakariakov] 
2009; Zimovets & Struminsky 2010; Nakariakov et al. 2010 J 
who claim a significant QPP detection, applied a standard Lomb- 
Scargle analysis to detrended solar flare light curves. Our inves- 
tigation of this method suggests the power in the low frequency 
range is being artificially suppressed, which can lead to mislead- 
ing values for the significances of features in the PSD. A peri- 
odogram analysis should always be performed on the raw and 
undetrended light curve as was done here. 
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We stress once more that not only solar flares but many astro- 
physical sources (X-ray binaries, Seyfert galaxies, GRBs) suf- 
fer from steep power spectra in the low-frequency range. Such 
power spectra make a periodogram analysis not trivial, the inter- 
pretation of a peak in a PSD more complex and an estimation 
of its significance an important issue. Again we emphasize that 
red-noise is an intrinsic source property. In other words, having 
shown that the variations in the corresponding solar flares are 
not quasi-periodic at the 3cr level does not mean that these vari- 
ations are not real. The corresponding flux changes are some- 
times dramatic, reaching a factor of a few within a few tens 
of seconds. Also, these variations occur in phase at different 
X-ray to gamma-ray energies, and other flares have been ob- 
serv ed to also occur in phase with microwave r adio emission 
(e.g. lNakariakov et al.ll2010al:lFoullon eTani2010h . 

While we think that it is not ne cessary to invoke oscilla- 
tory regimes of plasma inst abilities dNakariakov et al.l 1201 Obt 
iReznikova & Shibasakill201 1,) . we note that the physics of these 
variations is certainly interesting and worth further studies. 
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